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We have performed two kinds of non-linear fitting procedures to the single-particle energies in the 
sdgh major shell to obtain better shell model results. The low-lying energy eigenvalues of the light 
Sn isotopes with A = 103 — 110 in the sdg/i-shell are calculated in the framework of the nuclear 
shell model by using CD-Bonn two-body effective nucleon-nucleon interaction. The obtained energy 
eigenvalues are fitted to the corresponding experimental values by using two different non-linear 
fitting procedures, i.e., downhill simplex method and clonal selection method. The unknown single¬ 
particle energies of the states 2 si/ 2 , 1 ^ 3 / 2 , and 0 /iii /2 are used in the fitting methods to obtain better 
spectra of the 104 , 106 , 108,110 isotopes. We compare the energy spectra of the I04.i06,i08,ii0gj^ 
I03,i05,i07,i09gj^ isotopes with/without a nonlinear fit to the experimental results. 

PACS numbers: 21.60.Cs, 21.10.-k, 21.60.Fw 


I. INTRODUCTION 

The nuclear shell model assumes any nuclei as a many- 
body quantum mechanical system specified by the Hamil¬ 
tonian in the second quantization 

H = Y^ e^aiaa + X 

The solution of this Hamiltonian for a series of nuclei in 
the periodic table implies deriving systematically some 
nuclear properties such as the energy spectra, electro¬ 
magnetic transitions, and beta decay, etc [H, 0- 

The quantum many-body solution of the Hamiltonian 
in Eq. o for a series of nuclei is very complex that 
is normally divided in two sub-problems: the first one 
being the establishment of an effective interaction 
from the bare nucleon-nucleon (NN) potential in a va¬ 
lence space and of a set of valence space single-particle 
energies for sequences of nuclei. The second problem is 
the quantum many-body calculations of the Hamiltonian 
in Eq. O i- In this article, we present an improvement 
under the first problem. 

The one-body part of the Hamiltonian in Eq. O, i.e., 
the first term, describes the interaction of the valence 
nucleons with the closed core. The second term is the 
two-body part which describes the two-body effective in¬ 
teraction. In this study we have used the CD-Bonn two- 
body effective interactions in the sdgf/i-shell provided to 
us by M. Hjorth-Jensen et al. [3|. They calculate the 
charge-dependent Bonn (CD-Bonn) two-body effective 
interaction for shell model studies based on the free 
nucleon-nucleon interaction V However, the single¬ 

particle energies can not be obtained via their procedure. 
The latter are normally taken from the experimental data 
on the nuclei one particle away from the closed shell nu¬ 
clei. Unfortunately, in the sdgh major shell, the nuclei 
one particle away from the closed shell nucleus ^°°Sn have 


so far escaped by measurements. Therefore, we used in¬ 
stead the single-particle energies for neutrons in the sdgh 
major shell in our earlier studies around H = 100 mass 
region [mi, e.g., eid,/, = 0.00 MeV, = 0.08 MeV, 
£2*1/2 = 2-45 MeV, £1^3/2 = 2.55 MeV, and £0/111/2 = 3.20 
MeV. These energies are in reasonable agreement with 
similar shell model calculations in this region, see for ex¬ 
ample Ref. [l4jT6j. However, in the recent experimental 
work of Ref. [ll an excited state at 171.7 keV has been 
identified in ^°^Sn and interpreted as the 057/2 single¬ 
neutron state. 

The experiment in Ref. fixes the ambiguity about 
the single-particle energy of the O57/2 orbit. So, the un¬ 
certain single-particle energies in the sdgh major shell 
are left only on three of them, namely £2*1/21 £id3/2i 
and £0/111/2 • The ambiguity of these single-particle en¬ 
ergies establishes the need of this work. In this work, 
we have performed two kinds of non-linear fitting meth¬ 
ods, namely downhill simplex method and clonal selec¬ 
tion principle, to obtain the optimum single-particle en¬ 
ergies for the single-particle orbits of 2si/2, 1^3/2, and 
O/111/2. 

In order to obtain a better energy spectrum of any 
isotopes in sdgh major shell in the nuclear shell model 
framework we should improve both the two-body effec¬ 
tive interactions and the sin gle- particle energies given in 
Eq. ID). Wildenthal in Ref. [1^ provided a fit procedure 
in the sd-shell whereby be improved on an effective in¬ 
teraction by fitting the 3 single-particle energies and 63 
effective two-body matrix elements. The question being 
considered in that study is: What is the best possible 
shell model description of the nuclei in the sd-shell inde¬ 
pendently of what the Hamiltonian is? The same ambi¬ 
tious question was raised by Honma et al. [l^ for the 
p/-shell, where there are 4 single-particle energies and 
195 effective two-body matrix elements to fit. 

In the sdg/i-shell there are 5 single-particle energies 
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and 542 two-body matrix elements of effective interaction 
that specifies the Hamiltonian. A global fit of both 
the single-particle energies and these two-body matrix 
elements is certainly out of question. Many of the sdgh- 
shell nuclei can not be calculated yet, some of the nuclei 
that can be modeled require very large calculations and 
we have limited experimental data. This study restricts 
itself to fit single-particle energies. 

We have five single-particle states in the sdgh ma¬ 
jor shell. All theoretical and experimental data indicate 
that the lowest single-particle state is the 1 ^ 5/2 state and 
the energy of the 037/2 state are experimentally known 
saiMS- So we fix the 1 ^ 5/2 state as the lowest at the 
energy of 0.00 MeV and the O 37/2 state as the first ex¬ 
cited state at the energy of 0.172 MeV, and then perform 
a non-linear fit to the other three singlemarticle energies. 
We use the downhill simplex method and the clonal 
selection principle pH ] in three dimensions to fit these 
three unknown single-particle energies. The downhill 
simplex method was chosen over the others, i.e., steep¬ 
est decent annealing, monte carlo, because of the claim 
that it avoids trapping in local minima. The clonal selec¬ 
tion principle was chosen because it represents a global 
fit over the all cells. 


II. FITTING SINGLE PARTICLE ENERGIES IN 

3-D 

Because of the complexity and very large model space 
of the two-body effective interaction in the sdgh ma¬ 
jor shell, we only apply a fit process to the values of 
the single-particle energies. We adjust the values of the 
single-particle energies so as to fit the theoretical energy 
levels of the selected nuclei to the experimental energy 
levels. We outline some major steps in the present study 
and the fitting procedures below. 

For a set of known experimental energy data = 

Ji, J 2 ,...) with the angular momentum J of a specific nu¬ 
cleus with a mass A, we calculate corresponding shell- 
model energy eigenvalues . We minimize the 

quantity 

E (Eip-ELeoryf (2) 

A J=Ji,J2,... 

by varying the values of the single-particle energies. Since 
the fitting methods that we use, simplex and clonal se¬ 
lection, are the non-linear processes with respect to the 
parameters, it is solved in an iterative way with succes¬ 
sive variations of those parameters (single-particle ener¬ 
gies) followed by diagonalization of the Hamiltonian until 
convergence. 

We choose the light even-even Tin isotopes experimen¬ 
tally available. Since the fit process is a non-linear pro¬ 
cess, the selection of the heavier nuclei would require 
very long calculation times until convergence. Therefore 
the nuclei I04,i06,i08,ii0gj^ largest Tin isotopes for 


which we can fit the single-particle energies in some rea¬ 
sonable time. On the other hand, the more number of 
nuclei are included in the fit process, the more accurate 
fitted single-particle energies are obtained in the sdgh 
major shell. 

There are 5 single-particle energies in sdgh-shell model 
space. For simplicity we reduce the dimension of the 
problem to 3-D instead of 5-D by fixing the energies of 
the ld^f 2 state to be 0.000 MeV and the O 37/2 state to be 
0.172 MeV. The energy values of the other three single¬ 
particle states 2 si/ 2 , 1 ^ 3 / 2 , and 0 /iii /2 are varied during 
fitting process. The experimental energies used for the fit 
are limited to those of the first occurrences of low-lying 
states. The theoretical energies are also chosen to be the 
lowest in the theoretical spectra. 


A. The Downhill Simplex Method in 
Multidimensions 

The downhill simplex method is a geometrical hill 
climbing scheme to find the minimum of a function of 
more than one independent variable. The method re¬ 
quires only function evaluations. A simplex is a geomet¬ 
rical figure with V -|- 1 points (or vertices) that lives in 
a parameter space of dimension N and all interconnect¬ 
ing line segments and polygonal faces of vertices. In 2-D 
space a simplex is a triangle, in 3-D space it is a tetrahe¬ 
dron, and so on. Given the function value at each of the 
simplex’s vertices, the worst vertex is displaced by hav¬ 
ing the simplex undergo one of four possible “moves”, 
namely reflection, reflection and expansion, contraction, 
or multiple contraction (see more details in Ref. [^.[^b 

The downhill simplex method starts with V -|-1 points 
defining an initial simplex. If one of fV -|- 1 points is the 
initial starting point Pq, then the other N points are 

Pi = Pq + Xci (3) 

where ebs are N unit vectors, and A is a constant charac¬ 
terizing the problem’s length scale ( or A may be different 
on each direction). 

The downhill simplex method takes a series of moves, 
most moves just moving the point of the simplex where 
the function is largest (“highest point”) through the op¬ 
posite face of the simplex to a lower point. These moves 
are called reflections, and they are constructed to con¬ 
serve the volume of the simplex. When it can do so, the 
method expands the simplex in one or another direction 
to take larger moves. When it reaches a “valley floor” 
the simplex contracts itself in the transverse direction 
and tries oozing down the valley. If there is a situation 
where the simplex is trying to “pass through the eye of 
a needle” it contracts itself in all directions, pulling itself 
in around its lowest (best) point. The simplex under¬ 
goes successive such moves until no move can be found 
that leads to further improvement beyond some present 
tolerance. 


3 


The simplex method requires that one provides initial 
coordinates (x,y) for the simplex’s four vertices. Despite 
the simplex method’s pseudo-global abilities on a multi¬ 
dimensional and global problem, the choice of the initial 
location for the simplex often determines whether the 
global maximum or minimum is ultimately found. 


B. The Clonal Selection Principle 


Artificial Immune Systems (AIS) are computational 
methods inspired by biological principles of the natural 
immune system. AIS uses the ideas from immunology to 
develop an algorithm describing adaptive systems in vari¬ 
ous areas of science and engineering applications [H, [ 2 ^ . 
The clonal selection algorithm is based on the idea of AIS 
and focuses on a systematic view of the immune system. 
The algorithm is initially proposed to carry out machine¬ 
learning and pattern-recognition tasks and then is ap¬ 
plied to optimization problems 2M- optimization 

problems, this algorithm uses random instead of certain 
passing rules to avoid sticking on a local minima. 

Clonal selection algorithm uses two kinds of popula¬ 
tions of strings: a set of antigens (corresponding to a ob¬ 
jective function) and a set of antibodies (corresponding to 
variables of the objective function) within the maximum 
and minimum limits of variables. Antibody population 
is generated by random numbers scaled by limits of vari¬ 
ables. The affinity value of an antibody matches the value 
of purpose function calculated for a given antibody. The 
affinity population matrix is obtained by sorting affinity 
values. The amount of antibodies having highest affinity 
are selected from the entire population and a new set of 
population is formed as a sub-antibody array. 

A cloning ratio is determined by a percentage of each 
antibody array elements. Each of antibody array ele¬ 
ments is cloned by cloning ratio and inverting cloning 
factor determining the number of cloning. The affinity 
values of the new cloned population exposed to mutation 
are calculated and then a hypermutation process is per¬ 
formed by inverse cloning ratio. After the hypermutation 
process, a new number of antibody are determined and 
the cloned antibody array is obtained by sorting the cor¬ 
responding affinity values. An amount of antibodies with 
highest affinity are reselected and added to the antibody 
population set. Finally, an amount of antibodies with 
low affinity value within the antibody population set are 
replaced with the newly formed antibodies. At the end, 
the optimum values are determined. 


III. RESULTS AND DISCUSSION 

The initial values of the single-particle energies to be 
fitted are chosen to be £ 281/2 = 2.45 MeV, £ 1 ^ 3/2 = 2.55 
MeV, and £ 0 / 111/2 = 3.20 MeV. These energies are in rea¬ 
sonable agreement with similar shell model calculations 
in this region [ 1 ®. 
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FIG. 1: The single-particle energies before and after the fit 
processes. 


The calculation time for each iteration in both sim¬ 
plex and clonal selection fitting methods is very long 
time since it requires the shell-model calculations of the 
states J'^ = 0+, 2+, 4+, 6+, 8+, 10+ for the Tin isotopes 
I04,i06,i08,ii0gj^^ We have carried out the shell-model 
calculations in the High Performance Parallel Computer 
Cluster at Suleyman Demirel University, Turkey, consist¬ 
ing of 1 master and 30 slave nodes with each one of 3.2 
GHz Dual Xeon Processors. For a comparison, the re¬ 
quired time to calculate the energy of the state 4+ of 
^°®Sn and ^^°Sn are about 10 min. and 156 min., respec¬ 
tively. The largest Hamiltonian dimension among the 
states that we include in the fitting process is 156674 for 
the state 8+ of ^^°Sn. 

The fitted single-particle energies obtained by using 
clonal selection algorithm and simplex algorithm are 
shown in Figure [T] as compared to the ones before fit¬ 
ting. As mentioned above, the energies of the single¬ 
particle states ld 5/2 and 037/2 are fixed to the values of 
0.00 MeV and 0.17 MeV, respectively. Both fits give al¬ 
most same fitted energy of the single-particle state 2 si /2 
to be 1.56 MeV and 1.53 MeV by lowering almost 1 MeV 
from the initial value of 2.45 MeV. The energies of the 
























FIG. 2: Experimental and theoretical spectrum of and 

^°®Sn obtained with the single-particle energies before fit and 
the fitted ones using the CD-Bonn effective interaction. 



FIG. 3: Experimental and theoretical spectrum of ^°®Sn and 
^^°Sn obtained with the single-particle energies before fit and 
the fitted ones using the CD-Bonn effective interaction. 








FIG. 4: Experimental and theoretical spectrum of and 

^°®Sn obtained with the single-particle energies before fit and 
the fitted ones using the CD-Bonn effective interaction. 



FIG. 5: Experimental and theoretical spectrum of ^^’^Sn and 
^°®Sn obtained with the single-particle energies before fit and 
the fitted ones using the CD-Bonn effective interaction. 
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TABLE I: The a, root mean square (rms), values for the 2"^, 
4^, 6^, 8^, lO”*" states of Sn isotopes with A = 104-110 before 
and after fitting processes. 


+Sn 


^before 

•^clonal 

^simplex 


2 + 

0.16 

0.05 

0.02 


4+ 

0.05 

0.04 

0.06 

4°4Sn 

6 + 

0.03 

0.11 

0.11 


8 + 

0.01 

0.16 

0.16 


10 + 

0.09 

0.25 

0.22 


2 + 

0.16 

0.04 

0.02 


4+ 

0.14 

0.03 

0.00 


6 + 

0.06 

0.04 

0.03 


8 + 

0.11 

0.05 

0.05 


10 + 

0.15 

0.04 

0.00 


2 + 

0.25 

0.11 

0.03 


4+ 

0.27 

0.12 

0.07 


6 + 

0.12 

0.01 

0.01 


8 + 

0.32 

0.09 

0.08 


10 + 

0.16 

0.04 

0.00 


2 + 

0.32 

0.16 

0.06 


4+ 

0.34 

0.04 

0.06 

“°Sn 

6 + 

0.28 

0.13 

0.14 


8 + 

0.52 

0.06 

0.05 


10 + 

0.03 

0.28 

0.24 


other two single-particle states 1 ^ 3/2 and O/ 111/2 are very 
different each other after the fit processes. The simplex 
method changes the energy of the single-particle state 
1 ^ 3/2 -0.008 MeV as to be 2.47 from the initial value 
of 2.55 MeV while clonal selection method makes the 
change -1-0.65 MeV as to be 3.20 MeV. The major change 
happens to the energy of the single-particle state O/ 111 / 2 : 
with a change of -1-0.95 MeV and -1-1.6 MeV along with 
clonal selection and simplex methods, respectively. 

As we mentioned in Section II, the convergence of opti¬ 
mization algorithms are satisfied by minimizing the quan¬ 
tity, given in Eq. ©■ The value of tr^ before fit¬ 
ting process is 1.127. The values of cr^ after fitting pro¬ 
cesses are 0.208 after 76 iterations for simplex algorithm 
and 0.256 after 130 iterations for clonal selection algo¬ 
rithm. Both simplex and clonal selection algorithms quit 
the iterations by having unchange value of tr^, 0.208 and 
0.256 respectively. Therefore, the convergence of the fit¬ 
ted single-particle energies are obtained by lowering the 
value of cr^ from 1.127 to 0.208 and 0.256 for simplex and 
clonal selection methods, respectively. 

After having the fitted single-particle energies, we have 
carried out the shell-model calculations for the states O’*", 
2+, 4+, 6+, 8+, lO"*" of even-even ^°‘*d 06 .i 08 ,iiogj^ jgQ|.Qpgg 
to see how good spectrum obtained by using the new 
single-particle energies. In Table U we have given the 
root mean square (rms) deviations, tr, i.e., calculated by 
using Eq. [2 for the states 0+, 2+, 4+, 6+, 8+, 10+ of 
104406. 108 ,iiogj^ to see the goodness of the energy spectra 
before and after fitting processes. 

Figure [ 2 ] shows a comparison among the experimental, 
theoretical and fitted spectra of the i04Sn isotope. Both 


fits give very similar energies for the states 0+, 2+, 4+, 
6+, 8+, 10+ of i04Sn. Both clonal selection and simplex 
fits improve the energies of the low-spin states 2+ and 
4+ from 1.42 MeV to 1.31 MeV and 1.28 MeV, from 1.99 
MeV to 1.90 MeV and 1.88 MeV, respectively, compare 
to the experimental data. The energy of the states 6+, 
8+, 10+ is not improved through the both fits. 

Figure [2 also compares theoretical and fitted spectra 
of the i°®Sn isotope to the experimental spectra. All 
fitted states merge to the experimental ones. The states 
4+ and 10+ are exactly reproduced by simplex and the 
states 2+, 6+, 8+ are deviated slightly as to be 0.04, 0.04, 
0.05 MeV and 0.02, 0.03, 0.05 MeV by clonal selection 
and simplex methods, respectively. Overall, the simplex 
methods gives slightly better spectrum with respect to 
the clonal selection method for the 406Sn isotope. 

Figure [2 shows the theoretical (before fitting) and ex¬ 
perimental spectra of the ^°®Sn isotope as well as the 
fitted ones. From Table HI the energy deviations of the 
states to the experimental ones are 0.11, 0.03 MeV for 2+, 
0.12, 0.07 MeV for 4+, 0.01, 0.01 MeV for 6+, 0.09, 0.08 
MeV for 8+, and 0.04, 0.00 MeV for 10+ with clonal selec¬ 
tion and simplex methods, respectively. As we see from 
these rms deviations, all states merge gradually to the 
corresponding experimental ones. This mergence shows 
a definite improvement on the spectrum of the ^°®Sn iso¬ 
tope. 

Figure [3] also shows the theoretical (before fitting) and 
experimental spectra of the the ^i^Sn isotope as well as 
the fitted ones. All energy levels except the state 10+ 
merge to the experimental energies with the rms devia¬ 
tions of 0.16, 0.06 MeV for 2+, 0.04, 0.06 MeV for 4+, 
0.13, 0.14 MeV for 6+, 0.06, 0.05 MeV for 8+ of clonal 
selection and simplex methods, respectively. Both fit pro¬ 
cesses do not improve the energy level of the state 10+. 

If we look at the systematic behavior of the states com¬ 
pare to before and after fitting processes, the states 2+ 
and 4+ are definitely improved with both fitting proce¬ 
dures for all I04,i06,i08,ii0gj^ isotopes; the states 6+ and 
8+ are fitted well for the loe.ios.iiogj^ isotopes and not for 
4°4Sn ; the state 10+ is fitted well for i 06 ,i 08 gj^ isotopes 
and not for 404.iiogjj isotopes. 

To see the goodness of the obtained single-particle en¬ 
ergies we have also carried out the shell-model calcula¬ 
tions with the fitted single-particle energies for the low¬ 
est states 1^, I"*", I"*", I"*", of the light odd-even 
I03,i05,i07,i09gj^ isotopes. Figures m and [2 show such a 
spectrum comparing the calculated (before fitting) and 
available experimental values of the states § j § > ^ > 
of the i03.i05gj^ io7,io9gg isotopes, respec¬ 
tively, as well as the fitted ones. 

In Figure m the states and of 4^3 gg agj 7 + ^ 

I"*", of 433Sn are lowered well to get close to the 
experimental values. The states with higher energies of 
4°®Sn, e.g., I"*" and are fitted better with respect to 
the lower ones. 

In Figure[2the states I"*" of and ^^^Sn are fitted 
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to the values of 0.29 (0.27) and 0.08 (0.08) MeV with 
the rms deviations of 0.14 (0.12) MeV and 0.07 (0.07) 
MeV of clonal selection and simplex methods, respec- 
tively. The energy of | state is much better fitted to 
the corresponding experimental value by lowering about 
-0.20 MeV. As in the case of ^°®Sn, the states with higher 
energies of ^°^Sn and ^°®Sn are fitted better with respect 
to the lower ones. 

As we see from the Figure S] and [SJ the fitted spec¬ 
tra of odd-even I03,i05,i07,i09gj^ isotopes agree with the 
experimental data better than the ones of even-even 
104406, 108 ,iiogj^ isotopes. This is an evidence that the fit¬ 
ted single-particle energies give better spectrum for the 
other Sn isotopes which are not part of the fitting pro¬ 
cesses. 


IV. CONCLUSION 

The non-linear fitting procedures, namely clonal selec¬ 
tion and simplex methods, have been applied to optimize 
the single-particle energies in the sdgh major shell in 
the nuclear shell model calculations. The fitted single¬ 
particle energies are calculated as to be e 2 si /2 =1-53 
(1.56) MeV, £id 3/2 =2.47 (3.20) MeV, and =4-80 

(4.15) MeV with the simplex algorithm (the clonal selec¬ 
tion algorithm). The fitted energy of the single-particle 
state 2 si /2 are almost the same in both fit processes, 
but the other two fitted energies of the single-particle 
states are very different from each other. The optimiza¬ 
tion ideas of the simplex and clonal selection methods 
give this difference. The obtained fitted single-particle 
energies are still in the range of expectation. Indeed, 
there is a minimum and maximum limits of energy vari¬ 
ables within the clonal selection algorithm while there is 
no limits of variables within the simplex algorithm. 

For the point of shell model calculations, we have re¬ 


peated the shell model calculations of the states O"*", 2+, 
4 +, 6+, 8+, 10+ of the io4,io6,io8,iiogj^ isotopes with the 
single-particle energies obtained before and after fitting 
processes to get better agreement of shell model spectra 
with the experiment. We have obtained better energy 
spectra of the I04,i06,i08,ii0gj^ isotopes with the fitted 
single-particle energies. The fit processes improve the all 
states of the I04,i06,i08,ii0gj^ isotopes except the higher 
excited states of 404Sn and the state 10+ of 440Sn. So, 
the overall agreement of the fitted spectra is very good 
with experiment. 

We have also performed the shell model calculations 
of the states I"*", I"*", I"*", I"*", of the light odd-even 
103,105.107,i09gj^ isotopes with the single-particle energies 
obtained before and after fitting processes. We have ob¬ 
tained even better energy spectra for the I03,i05,i07,i09gj^ 
isotopes with respect to the I04,i06,i08,ii0gj^ isotopes 
which are parts of the fitting processes. 

The htted single-particle energies with the simplex 
method give the spectra slightly better corresponding to 
the experiment with respect to the ones with clonal se¬ 
lection method. The fitted single-particle energies can 
be used to do the shell model calculations for the other 
nuclei in the vicinity of H’lSn, especially for mid-heavy 
Tin (Sn) isotopes with A = 110 — 120 and light Antimony 
(Sb) isotopes. With this study, we present an improve¬ 
ment on the ambiguity of the single-particle energies to 
be used in future shell model calculations in the sdgh 
major shell. 
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